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Nuclei are propitious tools to investigate the role of the superfluidity in the compressibility of 
a Fermionic system. The centroid of the Giant Monopole Resonance (GMR) in Tin isotopes is 
predicted using a constrained Hartree-Fock Bogoliubov approach, ensuring a full self-consistent 
treatment. Superfluidity is found to favour the compressibitily of nuclei. Pairing correlations explain 
why doubly magic imclei such as ^"^^Pb are stiffer compared to open-shell nuclei. Fully self-consistent 
predictions of the GMR on an isotopic chain should be the way to microscopically extract both 
^\ ' the incompressibility and the density dependence of a given energy functional. The macroscopic 

, extraction of Kgym, the asymmetry incompressibility, is questioned. Investigations of the GMR in 

■ unstable nuclei are called for. Pairing gap dependence of the nuclear matter incompressibility should 

' also be investigated. 

. PACS numbers; 21.10.Re,24.30.Cz,21.60.Jz 



O 

(N 



The effect of the superfluidity on the compressibility of a Fermionic system remains an open question. Superfluidity 
initially referred to a system with a dramatic drop of its viscosity [l[ : it could be suspected that a super- fluid would 
be easier to compress than a normal fluid. This question has been investigated in Fermionic atoms traps, by studying 
the frequency of the compression mode with respect to the scattering length. Experimentally some increase of the 
■ frequency may be observed in the weak pairing regime 0, but this signal remains to be conflrmed. Theoretically 
I , , both microscopic and hydrodynamical investigations show no variation of the compression mode between the normal 
^ ■ and the superfluid phases [sj, but the analysis is complicated by the temperature change between the two phases. 
^ [ In nuclear physics, the study of the role of superfluidity in the compressibility can also be performed: the isoscalar 
^H^. Giant Monopole Resonance (GMR) is a compression mode, allowing to probe for related superfluid effects. Ideal 
tools are especiallyisotopic chains, where pairing effects are evolving from normal (doubly-magic) nuclei to superfluid 
• (open-shell) ones Moreover, the incompressibility of nuclear matter is a basic parameter in calculations describing 
^ ' neutron stars or supernovae, where superfluid effects are known to occur [5]. 

, Constraining the nuclear incompressibility modulus Kqo with experimental data on the Giant Monopole Resonance 
is a longstanding problem. The flrst relevant approaches have been performed following the work of Blaizot and 
[ Pearson [gl, 0] : only microscopic predictions of the GMR compared with the data could validate the Kqo value of the 
. I functional which was used. Using a fully self-consistent approach such as the Random Phase Approximation (RPA) 
^ ; dl or the constrained Hartree-Fock (CHF) method, the GMR data on ^ospb currentlyprovides Koo^^ 230 MeV in 
non-relativistic approaches [To] whereas K^o— 260 MeV is obtained for relativistic one The roots of this puzzle 
between the two approaches is still an open question, but it has been shown that the neutron-proton asymmetry 
dependence of the incompressibility, denoted as Kgym, plays also a role, as well as the density dependence of the 
functional [l^: a value of Kqo— 250 MeV could be extracted from ^°^Pb data, using a non-relativistic CHF method 
with a modifled density dependence of the functional. Therefore, it should be noted that K^o cannot be extracted 
from the measurement of the GMR in a single nucleus: several parts of the functionnal are tested simultaneously, 
Ci namely Kj-o, Visym and its density dependence. 

A previous study on the role of superfluidity on nuclear incompressibility has been performed, finding a negligible 
effect but the theoretical approach was not self-consistent. Indeed self-consistency is crucial since pairing effects 
are expected to be small in the GMR: this high energy mode is mainly built from particle-hole configurations located 
far from the Fermi level, where pairing do not play a major role. However giant resonances are known to be very 
collective [31 and pairing can still have a sizable effect on the GMR properties: around 10 % on the centroid, which 
is the level of accuracy of present analysis on the extraction of Koo [ifl \lM- This requires the advent of accurate 
microsco pic models in the pairing channel, such as fully self-consistent Quasiparticle Random Phase Approximation 
(QRPA) achieved only recently. Experimentally, the measurement of the GMR on an isotopic chain facilitates 
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the study of superfluidityon the GMR properties [17[, and the possibility to measure the GMR in unstable nuclei 
emphasizes this feature [18 ]. 

It is therefore necessary to go towards the measurement of the GMR on several nuclei, such as an isotopic chain. 
The overused method of precise GMR measurements in a single nucleus, such as ^°*Pb, may not be the relevant 
approach. Other nuclei have been used such as ^°Zr and ^^''Sm. Indeed when considering the available GMR data 
from which the K^o value has been extracted, ^°*Pb is stiffer than the Sn, Zr and Sm nuclei: Koo is about 20 MeV 
larger, both in non-relativistic and in relativistic approaches [13, [HI- The question may not be "why Tin are so soft 
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?"[T| but rather "why ^ospb is so stiff ?". 

Recently the GMR was raeasured on the stable Tin istopic chain (from ^^^Sn to ^^''Sn) Once again it has 

been noticed that it is not possible to describe the GMR both in Sn and in Pb with the same functional, Tin beeing 
softer than Pb [l^[l^. In the non relativistic case, fully self-consistent QRPA calculations on Sn isotopes lead to 
Kqo ~ 215 MeV. The relativistic DDME2 paremeterisation using QRPA describes well the GMR in the Sn isotopes 
20|, but predict a low value of the GMR in ^ospj-, compared to the experiment, as can be seen on Fig. 8 of Rcf. 

On the contrary, a recent relativistic functional describes well the ^"^^Pb GMR, but systematically overestimates 
the Sn GMR values of about 1 MeV [2l| . Finally, attempts have been performed in order to describe Sn GMR data 
with relativistic fuctionals having a lower incompressibility and hence different density dependence and K^sym value 
[121 . Once again the ^"^Pb and Tin GMR cannot be described at the same time: this puzzling situation is due to the 
higher value of Kqo extracted from ^°*Pb, compared to Tin, Sm and Zr nuclei. 

In Ref. [l^ it has been found that including pairing effects in the description of the GMR allows to explain part 
of the Sn softness : pairing may decrease the predicted centroid of the GMR of few hundreds of keV, located at ~ 16 
MeV. This is sufficient to change by about 10 MeV the extracted value of the incompressibility of nuclei Ka, defined 
as 0: 



J^GMR - \ Uj 

where m is the nucleon mass and (r^) is the ground-state mean square radius. 

In this work we follow up this idea and show that the consequences of superfluidity on nuclear incompressibility may 
solve the above mentioned puzzle. It should be noted that pairing is vanishing in the doubly magic ^ospj-, nucleus, 
unlikely the other nuclei. It is necessary to use a fully microsocopic method including an accurate pairing approach. 
In order to predict the GMR in a microscopic way we use the constrained HF method, extended to the full Bogoliubov 
pairing treatment (CHFB). The CIIF(B) method has the advantage to very precisely predict the centroid of the GMR 
using the m_i sumrule @, [l2|. The whole residual interaction (including spin-orbit and Coulomb terms) is taken into 
account and this method is by construction the best to predict the GMR centroid [l3|- Introducing the monopole 
operator as a constraint, the m_i value is obtained from the derivative of the mean value of this operator. The mi 
sumrule is extracted from the usual double commutator, using the Thouless theorem [2^. Finally the GMR centroid 
is given by ^GMR=\/rnrfrn~i. All details on the CHF method can be found in [ol [To|. 

The extension of the CH method to the CHFB case uses the HFB approach in coordinate space [13] • Skyrme 
functionals and a zero-range surface pairing interaction are used in this work. The magnitude of the pairing interaction 
is adjusted so to describe the trend of the neutron pairing gap in Tin isotopes. This interaction is known to describe 
a large variety of pairing effects in nuclei [25l| . 

Fig. 1 displays the GMR energy obtained from the Sn measurements (times A^/'^ to correct for the slow lowering 
of the GMR with the nuclear mass [i3|)- Microscopic CHFB predictions using two functionals are also shown: SLy4 
^ (Koo=230 MeV, which describes well the Pb GMR data), and SkM* [13 (Koo=215 MeV). Without pairing, the 
SLy4 interaction overestimates the Sn GMR data. Pairing effects (CHFB calculations) decrease the centroid of the 
GMR, getting closer to the data. This confirms the results of [l^, where a self-consistent HFB+QRPA approach was 
used to describe the Tin data. It should be noted that using a less general BCS approach for the pairing channel 
increases the GMR energy [l^, [11, HI] . This is due to the problematic treatment of high energy single-particle states 
in the BCS model. It is therefore important to use the full HFB approach to correctly describe pairing effects on the 
GMR, especially for the pairing residual interaction. Fig. 1 shows, as expected, that SkM* predictions with pairing 
effects are in better agreement with the Sn data. 

The striking feature of Fig. 1 is the peak of the GMR centroid, located at the doubly magic ^^^Sn nucleus, using 
the CHFB predictions. This indicates that pairing effects should be considered to describe the behaviour of nuclear 
incompressibility, and that vanishing of pairing make the nuclei stiffer to compress, confirming our previous statement 
on the stiffness of ^"^^Pb. The importance of pairing effect can be understood in a simple way: since the nuclear 
incompressibility is defined as the second derivative of the energy functional at saturation density [^, there is no 
obvious reason why the pairing terms of the functional would play no role in the nuclear incompressibility. Nuclear 
incompressibility is indeed very close from a residual interaction (as a second derivative of the energy functionals 
with respect to the density), and it is known that pairing effects are relevant in residual interactions [1^. This is 
straightforward in nuclear matter where Kco is expressed from the Fq Landau parameter [29]. However, on Fig. 1, 
the GMR centroid is shifted to lower energies for more neutron-rich nuclei than ^^^Sn, also because of the appearance 
of a soft L=0 mode, predicted in QRPA calculations beyond ^'^^Sn [soj . 

To further investigate the role of pairing on nuclear incompressibility. Fig. 2 displays (defined by Eq. (1)) 
with respect to the average pairing gap calculated using the HFB approach, from ^^^Sn to ^'^^Sn. A clear correlation 
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is observed: the more superfluid the nuclei, the lower the incompressibility. Hence it may be easier to compress 
superfluid nuclei. This may be the first evidence of the role of superfluidity on the compressibility of a Fermionic 
system. A possible interpretation is that Cooper pairs can modify bulk properties, as known from nuclear physics 
phenomenology 

The decrease of incompressibility in superfluid nuclei raises the question of a similar effect in infinite nuclear matter: 
for now, incompressibility is given independently from the pairing part of the functional. However, considering present 
results, equations of state used for neutron star and supernovae predictions should take into account pairing to provide 
their incompressibility value. The comparison with GMR data shows, as mentioned above, that the functional as a 
whole (including pairing effects) is probed. The question of the behaviour of Koo with respect to the pairing gap is 
raised: it seems clear from nuclear data that nuclei incompressibility K^a decreases with increasing pairing gap. This 
should be investigated in nuclear matter. 

It should be noted that recent attempts have been performed to extract the Kgy^ value, and its corresponding 
quantity in nuclei (K^): the incompressibility in the Sn isotopic chain has been studied using the macroscopic (liquid 
drop) formula of nuclei incompressibility derived by Blaizot @, [TtI . [sij . However the effect of pairing demonstrated 
above shows that the current macroscopic approach may not be well designed: on Fig. 2, pairing effects induce ^ 10 
MeV change on the nuclear incompressibility K^. Hence the macroscopic expression of nuclei incompressibility Ka 
should be extended to these terms, since the appropriate definition of nuclei incompressibility is the second derivative 
of the microscopic energy functional @. Presently, the microscopic approach is more relevant to extract Ksym'- the 
extraction of Koo, i^sym and the density dependence of the functional are related, as stated in Ref. The GMR 

data on isotopic chains should be used on this purpose. 

It is not possible to describe the GMR centroid of both ^°*Pb and other nuclei with the same functional, as stated 
above. The puzzle of the stiffness of ^"^^Pb may come from its doubly magic behaviour. In Fig. 1 there is a sharp peak 
at doubly magic ^'^^Sn, and it would be very interesting to measure the GMR in this unstable nucleus. It should be 
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FIG. 2: Nuclear incompressibilities Ka in "2-i32sn 

isotopes calculated with the CHFB method and the SkM* interaction, as 
a function of the pairing gap A predicted by the HFB calculation. 

noted that such experiments are now feasible A possible explanation of the ^osp]^ stiffness is that the experimental 
data of Fjgmr is especially increased in the case of doubly magic nuclei, as observed in ^ospj^ compared to the GMR 
data available in other nuclei (such as the Tin isotopic chain) . This difficulty to describe with a single functional both 
doubly magic and other nuclei has already been observed on the masses, namely the so-called "mutually enhancement 
magicity" (MEM) , described in [H, [s^l : functionals designed to describe masses of open-shell nuclei cannot predict 
the masses of doubly magic nuclei such as ^"^^Sn and ^°®Pb, which are systematically more bound that predicted. In 
order to consider MEM, it may be necessary to take into account quadrupole correlation effects due to the flatness of 
the potentials for open-shell nuclei [sj] . K^i being related to the second derivative of the energy with respect to the 
density, it would be useful to find a way to predict the GMR beyond QRPA by taking into account such quadrupole 
correlations. This may solve the current puzzle of the stiffness of ^"^Pb. It should be noted that it would also be 
relevant to measure the GMR on the Pb isotopic chain in order to provide a similar analysis than the one on the Sn 
nuclei. 

In conclusion, it is shown that superfluidity favours the compressibility of nuclei, using a fully microscopic CHFB 
approach on the Tin isotopic chain. This may be the first evidence of a sizable effect of superfludity on the compress- 
ibility of a Fermionic system. Pairing effects should be described using a full microscopic HFB treatment. Doubly 
magic nuclei exhibit a specific increase of the GMR energy, due to the collapse of pairing. ^°*Pb is therefore the 
"anomalous" data compared to the others. It is not possible to disentangle pairing interaction from the equation 
of state when providing the nuclear matter value of Kqo . Indeed the pairing gap dependence on the nuclear matter 
incompressibility should be investigated, since it is shown that incompressibility decreases with increasing pairing gap 
in nuclei. Additional theoretical investigations are called for in order to predict the GMR including the mutually 
enhancement magicity effect. The macroscopic extraction of Kgym may be ill-defined and should be extended to 
include pairing effects. Experimentally, measurements of the GMR in unstable nuclei should be performed in doubly 
magic ^"^^Sn, as well as extending the measurement on the Sn and Pb isotopic chains. 

The author thanks G. Colo, M. Grasso, J. Margueron, Nguyen Van Giai, P. Ring, H. Sagawa and M. Urban for 
fruitful discussions about the results of this work. 
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